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Abstract 

Flow and Diffusion Distributed Structures (FDS) are stationary spatially periodic 
patterns that can be observed in reaction-diffusion-advection systems. These struc- 
tures arise when the flow rate exceeds a certain bifurcation point provided that 
concentrations of interacting species at the inlet differ from steady state values and 
the concentrations are held constant. Normally, theoretical studies of these patterns 
are developed without concerning a noise. In this paper we consider FDS for a more 
realistic conditions and assume that the inlet concentrations are perturbed by a 
small noise. When the flow rate is small, the FDS is linearly sensitive to noise at 
the inlet. Even weak fluctuations destroy the stationary pattern and an oscillatory 
solution appears instead. For higher flow rates the instability becomes nonlinear: 
the pattern remains unaltered for a weak noise and undergoes the destruction when 
the noise amplitude passes a certain threshold. We represent a detailed description 
of these effects and examine two scenarios for the stabilization. 
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1 Introduction 



React ion- diffusion system are known to demonstrate a variety of instabilities 
that result in stationary patterns or oscillatory behavior. Flow and Diffusion 
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Distributed Structure (FDS) is relatively new type of patterns that can appear 
in a reaction diffusion system in presence of an open flow. The necessary 
condition for this pattern to appear is oscillatory instability and constant, 
non steady state concentrations of interaction species at the inlet. If these 
condition are satisfied, the FDS appears when the flow rate passes above a 
certain critical point. Patterns of this type were first described by Kuznetsov et 
al. |1] and later they were studied in more details by Andresen et al. [2] . Soon 
after discovery these patterns were observed in the experiments by Ksein and 
Menzinger Later the experimental results were summarized in a survey 
paper [5]. 

There is a bit confusing variety of terms concerning this type of pattern for- 
mation phenomena. Initially, Ksern and Menzinger suggested the term "Flow 
distributed oscillations" (FDO) |i3j. Later, the abbreviation FDS was suggested 
by Satnoianu [Ef7] . In this paper we adhere to the last version and refer to 
the structures under consideration as FDS. (See also the discussion of the 
terminology in |8]). 

Extensive theoretical studies of new structures are provided in piSfTfS] . In par- 
ticular, these papers explore the relationship between FDS and other known 
types of instabilities in reaction-diffusion-advection systems: Turing patterns, 
DIFI and Hopf instabilities. As shown by McGraw and Menzinger in [8], the 
FDS (they use the term FDO, see the remark above) is closely related to the 
Hopf instability and DIFI, while the Turing pattern has a different nature. 

As mentioned before, an FDS pattern emerges provided that a system demon- 
strates an oscillatory instability. Often a convective instability of the oscil- 
latory mode is declared as a necessary condition for FDS: arbitrary small 
constant perturbation to the steady state at the inlet can grow to an FDS. 
But as shown by Kuptsov [lOj, the FDS can also appear in an absolute insta- 
bility domain. In this case, however, the transition to FDS is rigid. If the inlet 
perturbation is small, the oscillatory solution dominates an FDS. However, 
the increasing of the perturbation gives an advantage to the FDS and it grows 
and suppresses the oscillations. 

An experimental study of an FDS in presence of a differential transport is 
provided by Miguez et al. [TTIT2] . They consider a remarkable variety of im- 
pacts on an FDS pattern. These are a parallel formation of two neighboring 
FDSs, an FDS in a system with monotonically increasing flow rate, an FDS in 
2D system with sinusoidally varying ^/-boundary and also interaction of two 
perpendicular FDSs. This analysis reveals a high degree of robustness and 
structural stability of this type of patterns. 

A study of forcing to FDS is also presented in a work by Kuptsov et al. in [13] . 
This paper provides a numerical analysis of an FDS that is perturbed by a 
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small particle dragged by a flow. Papers [Hll4fl5j are devoted to the study 
of an FDS in presence of periodic oscillations at the inlet. We shell discuss 
these papers below in the flnal section, because they are closely related to our 
analysis. 

Recently, Kuptsov at al. represented a theoretical analysis of FDSs in a 2D 
system with Poiseuille flow [16]. Two basic types of patterns were observed. 
The flrst one is y-dependant version of ID FDS. The pattern consists of curved 
stationary stripes that are transverse to the flow direction. Patterns of this 
type were also observed in the experiments by Miguez et al. [I2j. The second 
type of structures is speciflc to the 2D case: the structure consists of several 
stationary longitudinal stripes. If both types of basic structures are allowed, 
a more complicated pattern appears as a combination of transverse and lon- 
gitudinal stripes. 

All previous theoretical approaches have been developed for perfectly constant 
condition at the inlet. But fluctuations, that are always present in natural sys- 
tems, can influence the pattern formation. As observed in many experiments, 
see for example [5] and ^7j, FDS near the critical point can be suppressed by 
oscillatory solutions induced by the inlet fluctuations. 

The effect of an inlet noise on convectively unstable distributed systems was 
studied by many authors. We recall that in the convectively unstable system 
a growing perturbation drifts downflow but decays when observed from any 
flxed point. By contrast, absolute instability implies that the perturbation, 
being initially localized, grows at any flxed point in space [T8][T9|IT||20] (see also 
reviews therein). In view of the extreme sensitivity of the convectively unstable 
state to perturbations, Deissler [21j, Deissler and Farmer [2^, and Borckmans 
et al. [23j conclude that wave-patterns may be generated and sustained due to 
ampliflcation of fluctuations from an upstream noise source. More subtle than 
a straightforward ampliflcation of noise, an interplay between fluctuations and 
dynamics is suggested by Landa [21] who treat the onset of turbulence in flow 
systems as a kind of noise-induced transition. 

Kuznetsov [23] studied the effect of inlet noise on a spatially uniform reaction- 
diffusion-advection system. (This problem is similar to the ours, but we also 
apply a constant perturbation at the inlet to obtain an FDS pattern.) A noise- 
induced absolute instability was observed: the inlet noise increases the critical 
flow rate for transition from absolute to convective instability, and the new 
critical value is proportional to the square root of the noise amplitude, i.e., it 
obeys a power law. 

In our previous paper [26] we already discussed the stability of FDSs to an inlet 
noise. In particular, we observed that when the flow rate is close to the FDS 
critical point, the noise, regardless of its amplitude, destroys the structure, 
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i.e., the FDS is linearly unstable. For higher flow rates the system passes a 
new critical point, thus becoming insensitive to a weak noise. Beyond this 
point, noise destroys the structure only when the noise amplitude exceeds a 
threshold value. The threshold grows as a power of the flow rate. We shall 
refer to this case as nonlinear instability of an FDS. 

The purpose of the present paper is to extend our previous study. We demon- 
strate two scenarios of stabilization and develop a linear stability analysis of 
an FDS pattern. Sec. [2] reviews a variety of transitions to FDS. Sec. [3] repre- 
sents a qualitative picture of the destruction and stabilization of FDSs by a 
noise. In Sec. H] we summarize previous considerations and, employing Fourier 
analysis, introduce two scenarios of the stabilization. In Sec. this two sce- 
narios are confirmed by a linear stability analysis. Finally, Sec. |6] summarizes 
our results. 



2 Oscillatory instabilities and transition to FDS 

Before starting the main analysis let us introduce our model system and briefly 
review the oscillatory instabilities that can produce FDS patterns. 

We shall consider the well-known Lengyel- Epstein model for the CIMA (chlorite- 
iodide-malonic acid-starch) reaction [271l28ll6] . assuming that the reagents are 
carried by a flow in a one-dimensional reactor: 

Ut + (j)Ux — Uxx = a — u — Auv/ (l + u^^ , (la) 
vt + r(f)Vx - 6vxx =b(u- uv/ {l + u^)) ■ (lb) 

Here u and v denote the dynamical variables related to the concentrations of 
two chemical species, a and h are treated as control parameters, and (j) is the 
flow rate. The transport rates for u and v can be different: 5 introduces the 
differential diffusion, and r controls the differential flow. We shall study two 
cases: no differential flow at r = 1 and identical ratios for diffusion and flow 
at r = (5. The reactor is supposed to be semi-infinite. 

Note that we use a nonstandard form of Eqs ([1]). Normally, in the model equa- 
tions for CIMA reaction parameters 5 and h are proportional to each other. 
But in all our studies we always take a constant values for the parameters and 
analyze the varying of the flow rate or noise amplitude. So, at each step of ob- 
servations one can easily recalculate our parameters to the standard ones. But 
we believe that the specific numerical values are not very important because 
we do not study the system near bifurcation points. As we have examined, the 
represented results remain valid for a wide variety of the parameter values. 
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We shall find numerical solutions to Eqs. ([T]) via the semi-implicit Cranck- 
Nicholson method with steps of discretization Ax ?a 0.1 and At ^ 0.0125. 
Also, in all our numerical simulations an outlet condition {ux,Vx)x=x = is 
applied, where X is the length of the reactor. This type of the outlet conditions 
are normally used in numerical simulations of reaction-diffusion systems, see 
for example [28|IT] . These conditions are the best choice for a model of a semi- 
infinite, i.e., very long, system. In such a system the right boundary should not 
influence an overall picture. The vanish of and at the outlet effectively 
prevents an upstream propagation of any numerical artifacts. They, if appear 
at all, stay near the right boundary as a tiny area of oscillations. 

The system ([1]) admits the homogeneous steady state solution 

Us = a/5, t;s = l + aV5'. (2) 

In the absence of a flow the steady state can be destabilized either by Hopf or 
by Turing modes whose growth rates becomes positive at 

b = bn = (Sa^ - 125)/(5a) (3) 

and 

6 = 6t = (125 + 13a^ - 40^^10(25 + a'^))5) / {5a) , (4) 

respectively [2^30|28j . The Hopf instability gives rise to uniform temporal 
oscillations while the Turing instability produces a stationary space-periodic 
pattern. Note that for the considered system supercritical values of b lie below 
the bifurcation points. 

When both components of the system are transported by a flow with identical 
rates, i.e., r = 1, the Hopf instability remains unaltered while the Turing 
structure loses its stationarity. (Of course, this takes place only for an observer 
in a laboratory reference frame. The other one that moves with the flow still 
registers normal Turing stripes). The frequency of temporal oscillations of the 
appearing waves is proportional to the flow rate |31j . 

Differential flow with r ^ 1 results in the formation of new spatio-temporal 
structures due to the so-called Differential-Flow Instability (DIFI) [32] • The 
most interesting situation is observed outside the Hopf and Turing domains. 
The homogeneous steady state of the system remains stable until the flow rate 
reaches a critical point 0difi- Above this point, a DIFI mode become unstable 
and give rise to travelling waves. 

In presence of a flow an instability can be either absolute or convective. The ab- 
solute instability appears for small flow rates and manifests itself as a spreading 
of an initially localized perturbation over the whole space. For higher flow rates 
the instability becomes convective: the growing perturbation, being carried by 
the flow, decays in any fixed point. A critical point 0ca for the transition from 
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absolute to convective instability can be obtained from dispersion equation of 
a system as described in [T8l[T9l[T|20] (see also reviews therein). 

Figure d] represents solutions to Eqs. ([T]) that are induced by the instabilities 
described above. In the other words, this figure demonstrates what happens 
in our system without an FDS. For all panes of the figure the initial condi- 
tion at t = is the homogeneous steady state perturbed by small, spatially 
inhomogeneous fluctuations. The inlet is held at the steady state. 

Figure [11(a) illustrates a convective Hopf instability. We observe a temporally 
periodic structure with a wave number that is zero on average. Because the 
initial distribution of the reacting species is inhomogeneous, Hopf oscillations 
appears in different points with different initial phases. Zero wave number of 
a Hopf mode means that the effective interaction between neighboring points 
is very weak. As a result we observe in the figure a lot of local phase perturba- 
tions. The perturbations, however, are slowly smoothed out by the diffusion. 
These perturbations as well as the upstream edge of the structure travel down- 
flow. This confirms the convective nature of the observed instability. Note that 
this is dictated by a dispersion equation and does not depend on boundary 
or initial conditions. One can easily reproduce this convective moving of the 
spatio-temporal structure for different conditions at the system edges or for 
another inial perturbation. 

Figure [U^b) represents a spatio-temporal structure that appears due to an 
absolute Turing instability. The initial state for this figure is inhomogeneous, 
but, in contrast to the Hopf case, the structure forms due to the spatial in- 
teraction stimulated by a differential diffusion. This process subordinates the 
concentrations of the reacting species in different points and results in the 
perfect spatially periodic pattern. In a spatio-temporal diagram this should 
looks like a vertical stripes. But a flow transports the Turing structure. As a 
result we observe that the stripes are sloped exactly along the flow. (In the 
spatio-temporal diagram = 0.5, tmax = 200 and Xmax = 0^max = 100, so that 
the flow is parallel to the main diagonal). 

Finally, Fig. [T]^c) illustrates a DIFI structure. For the parameters values that 
are used for this figure, the system is stable without a flow, but different flow 
rates of the reactants produces an instability. Now both frequency and wave 
number are nonzero. Observe that they have opposite signs, and DIFI waves 
travel upstream. The DIFI takes place at (/>difi that is normally higher than 
the critical value for absolute-convective instability transition 0ca- In our case 
the DIFI is also convective, so that the structure is washed by a flow. 

FDS patterns appear in the system when temporal oscillations are spread in 
space by a flow. Every volume of an oscillating medium, being carried by the 
flow, appears with different phases in different spatial points and thereby pro- 
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Figure 1. Oscillatory instabilities in the system ([T]). Grey levels indicate values of u. 
For all panels a = 20. (a) Hopf instability at 6 = 6, 5 = 3, r = 1, (/> = 3. (b) Turing 
instability at 5 = 6, (5 = 15, r = 1, = 0.5. (c) DIFI at 6 = 12, 5 = r = 3, (/> = 1.5. 
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duces a spatially periodic structure. The flow rate is constant, so this structure 
is stationary if all portions have identical initial phases. For this to be fulfilled, 
a constant boundary perturbation should be applied at the inlet. This is a per- 
manent displacement of reagent concentrations from the homogeneous steady 
state. A wave number of FDS can be roughly estimated as cj/^, where uj is the 
frequency of the oscillations registered in a reference frame that moves with 
the flow (see discussion in [2|3|33|34] ). 

A Hopf mode is responsible for the formation of FDS when all components of 
the system have identical flow rates. In the presence of a differential flow this 
role is played by a DIFI mode. Turing oscillations can not produce an FDS 
because this solution is stationary in the co-moving reference frame, so that 
for a stationary observer every portion of a medium moves along the system 
with constant phase. This simple conclusion is in a good agreement with the 
analysis provided by McGraw and Menzinger [8]: the FDO/FDS patterns are 
shown to be closely related to Hopf instability and DIFI, while the Turing 
pattern formation mechanism has a quite different nature. 

To realize the formation of FDS, it is enough to perturb only one of the 
components and keep steady state value for the other: 

m(x = 0,t) = Ms + Mbnd, w(a; = 0,t) = ws, (5) 

where Wbnd is a constant perturbation of the boundary value. As before, Ms 
and f s represent the homogeneous steady state. When the flow rate is small, 
the perturbation decays in space. A non-decaying FDS solution appears above 
a critical point p!|2|6] . 



ma%{5 + rf - {2,a^5 + bah - 125^) V 
(25 + a2)(5 + r)(3a2r - bab - 125r)r ' ^ 

This critical value diverges at 6 = rh-^. Without a differential flow, i.e., at 
r = 1, the condition of divergence coincides with the condition for the Hopf 
instability, while at r > 1 an FDS appears even if the flow-less system is stable. 

li (j) > <^FDS and </> > (/)ca, the transition to FDS takes place when an oscillatory 
mode is convectively unstable. Because in this case growing oscillations travel 
with the flow, any small inlet perturbation can freely grow to the FDS pattern. 
This case is referred to as a soft transition to FDS pUfTU] . Otherwise, if (/>fds < 
(j) < (pea., the oscillatory mode is absolutely unstable and tends to occupy the 
whole space. Now the FDS can grow only when the inlet perturbation is high, 
because it must overcome the competition from the oscillatory solution. This 
is a rigid transition to FDS [20p0] . In this paper we analyze FDS for > 0ca; 
i.e., we are interested in the case of the soft transition. 
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3 FDS in presence of inlet noise: a qualitative picture 

We introduce a noise via the following boundary condition at the inlet: 

U{x = 0,t)=Us+ Mbnd(l + eW), V{X = 0,t)= Vs. (7) 

Here, Ms and vs denote the homogeneous steady state and Wbnd is the constant 
boundary perturbation, ^(t) describes a noise with the uniform distribution 
in the interval [—1,1], and 6 controls an amplitude of the noise relating to 
Mbnd- The particular distribution of a noise is not very important. But, as we 
explain below, it should be sufficiently wide to contain required harmonics. 

In our previous work [2B] we have studied the stability of FDS patterns in the 
Hopf oscillatory domain. In this paper we consider the following four cases. 

(1) FDS appears in presence of the Hopf instability: r = 1 and < b < b^. 

(2) FDS is observed in presence of a differential flow, r > 1. The system is 
stable without a flow, i.e., & > 6h and b > 6t- Oscillations appear due to 
DIFI when the flow rate is above a critical point 0difi- 

(3) The transition to FDS takes place when both Turing and Hopf modes are 
unstable: r = 1 and 6 <^ 6h < ^t- The Turing mode has a higher growth 
rate, but the Hopf mode is also signiflcant. 

(4) FDS appears again in the presence of Hopf and Turing modes, but the 
system is close to the Hopf point: b < bn. Growth rate of the Hopf mode 
is small, so that this mode is much weaker than the Turing one. 

The analysis developed in [20] shows that in presence of an inlet noise the FDS 
is destroyed at small flow rates and appears only when the flow rate passes a 
certain point of stabilization. The stabilization flow rate is found to depend 
on the noise amplitude as a power law. For a flxed flow rate the FDS is stable 
if the noise amplitude is small and the pattern is destabilized and destroyed 
when the amplitude passes a certain threshold. In the present paper we refer 
to this phenomenon as a nonlinear instability of FDS. The picture differs for 
flow rates a bit above the FDS critical point 0fds- The threshold noise level 
absent and FDSs are linearly unstable so that any small inlet fluctuations 
destroy the pattern. We already reported this case in [26j. In Sec. [5] a detailed 
investigation is provided. 

3.1 Case 1: Hopf instability 

Figure [2] demonstrates an interaction between an FDS and inlet noise when 
only a Hopf mode is unstable. The constant part of the inlet perturbation 
is small, so that the FDS in this flgure has a boundary layer. A structure 
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Figure 2. FDS in the presence of inlet noise when a Hopf mode is unstable. In panel 
(a) the noise destroys the FDS pattern, while in panel (c) (see below) the FDS is 
stabilized at a higher flow rate. Flow rate grows from (a) ^ = 6 to (b) (p = 6.5 
and (c) = 7. a = 20, 6 = 6, (5 = 3, r = 1, Ubnd = 0.1, 6 = 0.05 (6h = 10.75, 
6t ^ 3.289). 

without a boundary layer is stabilized at smaller flow rate and has a number 
of peculiarities that we shall discuss below. 

In Fig. [2](a) the FDS pattern disappears after several periods and a new os- 
cillatory structure forms instead. This effect can be explained by a selective 
amplification of Fourier components of the noise. The maximum of this ampli- 
fication corresponds to a Hopf mode [26]. This results in a resonant excitation 
of this mode by the noise so that a growing oscillatory solution destroys the 
FDS. A corresponding Hopf structure without a noise is shown in Fig. [I](a). 
Note a high degree of similarity between noise and noiseless solutions. 

If an FDS has a boundary layer, i.e., Ubnd is small, the amplification basically 
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Figure 2. {continued) 

occurs inside this layer. A perturbation to the pattern always manifests itself 
at the end of the boundary layer regardless of the flow rate, compare Figs. [2]^a) 
and (b). If Wbnd is high and the boundary layer is absent, the noise also destroys 
an FDS pattern. But an effective amplification ratio of a fully developed FDS 
pattern is much lower, so that the length of an unperturbed FDS area can be 
very long. (We do not provide an illustration because, except for the boundary 
layer, the spatio-temporal diagram looks as that in Fig. El^a)). 

Figure[2]^b) illustrates an intermediate situation for a higher flow rate. Now the 
influence of the noise is weaker and the FDS is not totally destroyed. The over- 
all structure is a composition of vertical fragments of FDS and sloped stripes 
of an oscillatory solution. Observe a remarkable regularity of this composition. 
This interesting structure appears because both the stationary stripes and the 
oscillatory fragments are different manifestations of a Hopf solution. (Recall 
that FDS is, in fact, an oscillations that are spread in space due to a flow and 
diffusion.) This explains why the fragment are joined so well into continuous 
step-like stripes. 



This type of structures is a typical response of an FDS to a perturbation. 
Kaern and Menzinger |1] first observed them in a system with oscillatory in- 
let boundary The structures was referred to as pulsating waves. A more 
complicated zigzag patterns are reported by Taylor et al. [17]. Also, similar 
structures was described by Kuptsov et al. [13] for an FDS perturbed by a 
moving particle. In addition, coexisting stationary and travelling waves have 
been observed in a system without noise [7]. 

FDS without a boundary layer can also demonstrate spatio-temporal struc- 
tures composed of stationary and oscillatory fragments. An example is shown 
in Fig. [31 As mentioned above, for FDS without a boundary layer the ampli- 
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Figure 3. Spatio-temporal intermittency for an FDS without a boundary layer. The 
parameters are as in Fig. [2] except cf) = 6.1 and Ubnd = 2. Unlike Fig. [2|^b), the 
structure appears far from the inlet. 

fication area can be very extended. In Fig. [3] the FDS at a; < 1800 looks like a 
noiseless pattern. However, beyond this point noise demonstrates itself. Note 
that the flow rate in Fig. [3] is less then in Fig. [2]^b). In the other words, an 
FDS without a boundary layer is less sensible to noise and becomes stable at 
a smaller value of the flow rate. 

Further increase of the flow rate results in the growth of the stationary frag- 
ments and shrinking of the oscillating ones. Finally, all oscillating fragments 
disappear and we observe a pattern as in Fig. WLs). The structure here looks 
like a perfect FDS without a noise. 

Figure H] provides a more accurate veriflcation of the stabilization. This flgure 
shows spatial distributions of variances of temporal oscillations. The distri- 
butions are calculated for the diagrams in Fig. [2l Labels (a), (b) and (c) on 
the curves correspond to the panels of Fig. [2l All curves grow exponentially 
within the boundary layer of the FDS and nearly coincide there. In the other 
words, the boundary layer amplifles noise even when the corresponding fully 
developed FDS is stable. Behind the boundary layer the curves behave dif- 
ferent. Curve (a) continues to grow and reaches a saturation very fast. This 
corresponds to the oscillations observed in the right part of Fig. [2]^a). The 
intermediate curve (b) also grows, but does not tend to saturation. It varies 
slowly and irregularly near a relatively high value. This occurs because time 
series in this case are composed of fragments of two solutions, see Fig. |2]^b). 
They are irregularly arranged and their sizes are distributed within a wide 
range. Curve (c) conflrms the stabilization. The variance decays, therefore, 
the FDS remains stable even far downflow. 
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Figure 4. Spatial distributions of the variances of temporal oscillations before, curves 
(a) and (b), and after the stabilization of FDS, curve (c). The curves correspond to 
the panels in Fig. [2] with respective letters. 
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Figure 5. Intermittent time series of u recorded in a fixed spatial point for the case 
represented in Fig. [2l^b). Laminar phases correspond to fragments of FDS stripes. 

The combination of stationary and oscillatory solutions that is demonstrated 
in Figs.[2](b) and[3]can be treated as a spatio-temporal intermittency. A similar 
compound structure has been analyzed by Kuptsov et al. [13]. Figure [5] demon- 
strates an example of an intermittent time series. Variable oiu is observed from 
a fixed point in space: laminar phases, that correspond to fragments of FDS, 
are interrupted by bursts of oscillations. 

A characteristic feature of intermittency is a power-law divergence of the mean 
length of laminar phases as a bifurcation parameter approaches a critical point. 
In our case the flow rate controls the mean length of FDS fragments and the 
divergence corresponds to the stabilization of FDS. Let us employ the following 
algorithm: Given the flow rate, we compute a pure FDS pattern without a 
noise. Then we take a point at some distance from the inlet where the FDS 
is definitely fully developed (for all computations we take the same point). 
Starting from this position, within the period of the FDS we seek another 
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Figure 6. Mean lengthes of FDS fragments vs. flow rate measured in a fixed spatial 
point. The parameters are as in Fig. [2j The dots, computed numerically, are plotted 
in double logarithmic scale. In the right part the numerical data are approximated 
well by a straight line that indicates a power-law dependence. The exponent, equal 
to the slope of the line, is j ^ 73. 

point where u is in the middle between its maximum and minimum values. 
The value of u in this point is treated as a laminar state. Now a noise is 
switched on and after a transient period the mean time of staying near the 
laminar state is registered. Repeating this procedure for different flow rates, 
we obtain the dependence shown in Fig. El When the flow rate is small, u 
oscillates, as shown in Fig. Et^a), and passes the laminar state very fast. In 
Fig. [6] the corresponding mean values are very small. Intermittency appears 
for > 6.2, and the mean time spent near the laminar state starts to grow with 
the flow rate. The numerical data in this part of the flgure are approximated 
well by a straight line in double logarithmic scale. This implies a power-law 
dependence < L >= {(j) — (pc)"^. Note that the exponent 7 ~ 73 is extremely 
high. 



3.2 Case 2: DIFI 

This case, illustrated in Fig. [3, is similar to the previous one. In Fig. El^a) 
we observe how a noise excites an unstable oscillatory mode that grows to an 
oscillatory solution. But now a DIFI mode is responsible for the formation of 
FDS. Hence, it is the DIFI solution that appears behind the destructed FDS. 
Observe the high regularity of the appearing structure. In Fig. [7](b) a higher 
flow rate results in the formation of spatio-temporal intermittency similar to 
the structure in Fig. [2](b). A perfect correspondence between stationary and 
oscillatory fragments is also observed. The explanation is as in the previous 
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Figure 7. FDS with noise when a DIFI mode is unstable: (a) cp = 2.3, and (b) 
= 2.47. a = 20, 6 = 12, 5 = r = 3, Ubnd = 0.1, 9 = 1 {bn = 10.75, br ~ 3.289, 
0DIFI ~ 1-47 and (ppj^s ~ 2.27). 

case: On the one hand, the DIFI mode produces the oscillatory solution and, 
on the other hand, it is responsible for the formation of FDS. Further increase 
of the flow rate results in the stabilization of the FDS. Note that the noise 
amplitude is very high: 6 = 1. This is even higher then the constant part 
Uhnd = 0.1 of the inlet perturbation. But this can not prevent the stabilization. 
The pattern looks as in Fig. M^c) and is not shown here. 

As one can see from Fig. [TJ the boundary layer of the FDS manifests itself in 
a similar way as in the case shown in Fig. m The boundary layer effectively 
amplifies the noise, so that the oscillatory solution appears near its right edge. 
A structure without a boundary layer can also amplify noise, but the area of 
amplification is very extended and grows when the flow rate becomes higher. 
In particular, one can observe spatio-temporal intermittency only far from the 
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Figure 8. FDS with noise in the presence of Turing and Hopf modes with significant 
growth rates. For all panels a = 20, 6 = 5, (5 = 20, r = 1, Ubnd = 0-5, = 0.1 
(6h = 10.75, 6t ~ 21.928). Panels correspond to different flow rates: a) = 10, b) 

(j) = 12. 

inlet. As in the previous case, FDS without a boundary layer is stabilized at 
a smaller flow rate in comparison to the structure that has a boundary layer. 
The described behavior of FDS with and without a boundary layer is typical. 
We have also analyzed the effect of boundary layers for the two remaining 
cases and obtained similar results. 



3.3 Case 3: Turing and Hopf modes with significant growth rates 



Now we assume that 6 <^ 6h < ^t, so the system displays both Hopf and Tur- 
ing significant unstable modes and the Turing mode has a higher growth rate. 
As follows from direct numerical simulations, in a corresponding unbounded 
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Figure 8. (continued) The flow rates in the panels are c) (j) = 14, d) (p = 15.5. 
noiseless system the Turing solution grows faster and suppresses the Hopf one. 

Let us begin from small flow rates. When the FDS is affected by a noise, the 
Turing mode suppresses both the FDS and Hopf solutions and we observe 
almost perfect Turing structure behind a narrow FDS FigEKa). This 

figure is qualitatively similar to Figs. [2]^a) and[7](a) where unstable oscillatory 
modes also grow to form spatio-temporal structures in place of a destroyed 
FDS. 

Contrary to the previous cases, the higher flow rate in Fig. [8](b) does not 
result in the reestablishment of FDS. It gives an advantage to the Hopf mode 
instead. This one manifests itself as localized travelling areas of oscillations 
embedded into the Turing structure. 

Further increase of the flow rate makes the Hopf mode stronger, see Fig. [8]^c). 
The Turing waves appear now as localized perturbations to the Hopf structure. 
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In the other words, we have another spatio-temporal structure produced by 
the strongest oscillatory mode that suppresses the FDS. 

If the flow rate continues to grow, the Hopf solution starts to give up its 
area to the FDS. This produces a spatio-temporal intermittency. Fig. M^d). 
Note that a Turing structures disappear in a different manner. In this flgure 
there are rare localized Turing waves that are swept by the flow and decay in 
space. Sufliciently far from the inlet a combination of oscillatory and stationary 
solutions appears in the same manner as in the previous cases, see Figs. [21(b) 
and[7l^b). As the flow rate increases even further a stabilization of the FDS 
takes place. We do not provide a corresponding figure because the pattern 
looks like in Fig. [2](c). 

3.4 Case 4- Turing instability in presence of a weak Hopf mode 

This case differs from three previous. Though a Hopf mode is unstable and 
give rise to the formation of FDS, there is another unstable mode, namely 
the Turing one, that has much higher growth rate. When the flow rate is 
small, the spatio-temporal structure is as in the previous case. Figure M^a) 
shows how a noise stimulates the growth of a Turing solution in place of 
the FDS while the Hopf mode does not manifest itself at all (compare with 
Fig. [8](a)). But now even when the flow rate becomes higher a Hopf solution 
does not emerge. The left end of the Turing area becomes indented, and we 
observe an intermittent structure. Fig. [9](b). Unlike three previous structures, 
oscillatory and stationary fragments here have different origins. The stationary 
FDS stripes appear due to the Hopf mode while the oscillatory segments are 
produced by the Turing mode. Hence, fragments of different solutions are not 
joined into continuous, smooth zig-zag structures. They overlap and partially 
suppress each other instead. When the flow rate continues to grow, the FDS 
segments are enlarged while the Turing areas becomes more and more rare. 
Finally, the FDS becomes stable and occupies the whole space. The resulting 
pattern looks like in Fig. [2]^c). 



4 Two scenarios of stabilization 

Four cases considered above demonstrate two different scenarios of stabiliza- 
tion. The flrst one involves cases 1, 2 and 3, while the second is represented 
by the case 4. A clearer evidence of this can be provided with the Fourier 
analysis: We are going to demonstrate the dependance of Fourier spectra of 
time oscillations on the flow rate. This can be done with the following proce- 
dure. A point is flxed near the outlet and a series of temporal oscillations is 
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Figure 9. FDS with noise when a Turing mode is strong while a Hopf mode is 
weak. Observe the formation of Turing structure in place of the FDS in panel (a), 
(j) = 22, and spatio-temporal intermittency in panel (b), (j) = 23.5. The parameters 
are a = 25, 6 = 10 (6h = 14, 6t = 37.549), 5 = 25, r = 1, Ubnd = 0.5, 9 = 10"^ 



recorded in this point. Then the mean value of the series is subtracted and 
Fourier spectrum is computed using the FFT algorithm. This is repeated for 
different flow rates and the resulting Fourier amphtudes are plotted via grey 
scales in Fig. [TOl The panels (a), (b), (c) and (d) in this figure correspond to 
Figs. [21 [71 [8] and [9], respectively. 

As we already discussed in our previous paper [20], the destruction of FDS 
takes place because this pattern can selectively amplify a noise and the peak 
amplification corresponds to the unstable oscillatory mode. This is confirmed 
in Fig. [TOl We see here that Fourier components are gathered around an unsta- 
ble modes of the system. The resonant interaction between selected harmonics 
and internal unstable modes results in the destruction of FDS. 
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Figure 10. Fourier spectra of temporal oscillations vs. flow rate at a fixed point near 
the outlet. Grey levels indicate values of Fourier amplitudes. Panels (a), (b), (c) and 
(d) (see below) correspond to Figs. [21 [3 [8] and [9l respectively. Dashed lines show 
Hopf (H), Turing (T) and DIFI (D) frequencies. Light areas on the right parts of 
the figures correspond to the stable FDS. 
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Figures [TO](a) and [TOT b) demonstrate that the stabihzation in the Hopf and 
DIFI cases occurs due to the detuning of the selected harmonics from the 
resonance. Bands of selected frequencies are shifted down. Though the noise is 
still amplified, but the selected frequencies differ from the resonant and, hence, 
can not produce intensive excitation of the oscillations. The amplitudes of the 
outlet oscillations are reduced. We observe in the figure that the shaded areas 
becomes lighter as the flow rate grows. This corresponds to an intermittency 
in the spatio-temporal diagrams, see Figs. [21(b) and El^b). The shift of the 
frequencies is reflected in these figures by the partial straightening of the 
oscillatory segments that occur due to the merging with vertical FDS stripes 
(compare, for example, the Hopf structure in Fig. [2]^a) and the oscillatory 
areas in Fig. [2](b)). 

Figure [TUl(c) illustrates the same scenario of the stabilization. This takes place 
via the detuning of selected frequencies from the resonance with the Hopf 
mode. But here another band of Fourier components appears around the Tur- 
ing mode. When the flow rate is small, noise is not amplified near the weak 
Hopf mode. The resonance with the Turing mode is responsible for the for- 
mation of the oscillatory solution in Fig. [H](a). For a higher flow rate both 
modes are excited by the noise harmonics and we observe a spatio-temporal 
structure as in Fig. [8](b). Unlike the Hopf mode, the selective amplification 
around the Turing frequency does not depend directly on the flow rate and 
the Fourier components remain gathered around the Turing mode. But the 
coefficient of amplification vanishes as the fiow rate grows: the shaded area 
in the figure becomes lighter. The corresponding spatio-temporal structure is 
shown in Fig. [8](c). Finally, the Turing band disappears and further stabiliza- 
tion proceeds as in previous panels (a) and (b) of Fig. [TOl The corresponding 
spatio-temporal diagram in Fig. Et^d) is also similar to the Hopf and DIFI 
cases. Figs. |2](b) and[7](b), respectively. 

Fig. [TOT d) illustrates the second scenario of the stabilization. The selective 
amplification takes place only around the Turing mode, while the Hopf mode 
does not manifest itself except the small area at ~ 24. The gathering of 
the Fourier components near the Turing mode is similar to the case 3, see 
Fig. [TUT c). But now this is the main mechanism responsible for the stabiliza- 
tion. Observe that the shaded area in Fig. [TOlfd) disappears very sharply. This 
is explained by an intermittent nature of analyzed time series, see Fig. IHt^b). 
When the fiow rate approaches the stabilization point, intervals between Tur- 
ing bursts grow and sooner or later we get an interval that is longer than 
the observation time. In this point the characteristic Fourier spectrum sud- 
denly disappears. Though in the previous case 3 the Turing oscillations also 
demonstrate a kind of intermittency with the Hopf solution, see Fig. Mjo), 
the corresponding Fourier components in Fig. [TOT c) disappear smoothly. It 
indicates that there the intervals between Turing areas do not diverge. 
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5 Linear stability analysis 



In our previous paper [26] we have found that near the critical point 0fds an 
FDS pattern is linearly unstable and any small perturbation applied at the 
inlet can destroy it. In this section we study this effect and develop a linear 
stability analysis of the FDS. 

Let us first recall the discussions of Fig. HI This figure displays spatial distribu- 
tions of the variances of temporal oscillations of u. The variances, both before 
and after the stabilization point, grow exponentially within the boundary layer 
of FDS. The stabilization reveals itself only for a fully developed FDS while 
the boundary layer always amplifies noise. We have observed this for many 
parameter values and believe that this is typical for FDSs. Thus, we can ne- 
glect the boundary layer and consider a fully developed periodic FDS solution 
upDsi^ + -P) = '^FDsi^) and fFDs(^ + -P) = vpDsix), where P is the period. 
To analyze the stability of this solution we introduce a small perturbation, 
u = ufds{x) + [/(x)e~*'^* and v = wfds(^) + V(x)e~*'^*, where the frequency 
u plays a role of an additional parameter, u is real because the perturbation 
appears due to an external forcing that has a stationary amplitude. After the 
substitution into Eqs. ([1]), for small U{x) and V{x) we obtain the following 
linear equations with periodic coefficients: 

U" - (j)U' = {1 ~ iuj + AFi{x))U + 4^2 (8a) 
6V" ~ r(f)V' = b{Fi{x) - 1)U - {luj - bF2{x))V, (8b) 

where Fi = vpBsi^ - ^FDs)/i'^ + "fds)^ and F2 = mfds/(1 + u^j^s)- These 
second order equations have complex coefficients. After straightforward trans- 
formations they can be re-written as eight real equations of the first order. 

Similar approach in a context of FDO/FDS patterns was developed by Mc- 
Graw and Menzinger [S] who studied small perturbations to the steady state 
and considered linear modes having real frequencies. This type of modes cor- 
responds to a stationary forcing at the inlet. The essential difference of our 
study is the consideration of small perturbations to a fully developed FDS pat- 
tern. The linearized equations ([8]) in this case have periodic coefficients and 
we should apply the Floquet theorem to analyze its stability properties. 

As follows from the Floquet theorem ^35j the system has eight eigensolutions 
//(x) and eight corresponding eigenvalues A for which 

z7(x + P) = Az7(x). (9) 

There are four couples of A, either real and identical or complex conjugated. 
Their absolute values can be labelled as |A|i < |A|2 < |A|3 < |A|4. These values 
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depend on the parameters of the system and on the forcing frequency uj. For 
> 0FDS the two highest values always correspond to diverging solutions of 
Eqs ([HD, I Ala > 1 and |A|4 > 1, while the others always represents a stable 
solution, |A|i < 1. The stability of FDS is characterized by |A|2. FDS is stable 
if |A|2 < 1 for any lo while it is unstable when |A|2 > 1 for certain values of 
LO. It is more convenient to consider a characteristic exponent that does not 
depend on the period of FDS, 

Pi = (ln|A|i)/P. (10) 

So, in a critical point of the linear stabilization the global maximum of p2{}^^ 
passes zero. 

We compute the eigenvalue |A|2 using the Floquet theorem. The idea of this 
approach can be found, for example, in [55] . Let A be a linear evolution opera- 
tor advancing a solution over the period, Av[x) = v{x + P). Each eigenvector 
z7 can be decomposed over the unit basis as z7 = Yh=i eiQ, where Cj are the 
unit vectors and q are the coefficients of the decomposition. Taking this into 
account, we obtain from Eq. ([9]): 

8 8 
5] (Ciici) = A 5] (QCi) . (11) 
i=l i=l 

These equations have a nontrivial solution for the coefficients q if A is an 
eigenvalue of the matrix with columns Acj. To obtain this matrix, we find 
numerically a noiseless FDS solution and extract a full period. Then, we inter- 
polate it with polynomials and substitute the resulting functions into Eqs. ([H]). 
Starting from the unit vectors Cj (z = 1 ... 8), we find eight solutions to these 
equations over the period and obtain the sought matrix. Finally, we compute 
the eigenvalues and, after sorting by their absolute values, take the second 
one. 

A numerical analysis shows that on the FDS onset = 0fds the number of 
maxima of p2('^) coincides with the number of linear modes, stable or unstable, 
of Eqs. ([1]). There is at least one couple of symmetric maxima corresponding to 
Hopf or DIFI modes that are responsible for the formation of FDS. Without a 
differential flow, the system can also display a Turing mode. In this case P2('^) 
acquires another couple of maxima. For example in Fig. [11] the flrst maxima 
from the origin correspond to the Hopf mode and the second pair correspond 
to the Turing mode. All the maxima are positive, hence the FDS is unstable. 

The system ([S]) can be considered as a linear amplifler of an inlet forcing where 
p2{oj) is an ampliflcation ratio. As follows from Fig. [11] the ampliflcation is se- 
lective. When a wide-band noise passes thorough the system, noise harmonics 
grow if corresponding ^2(1^) is positive. All other harmonics decay. As a re- 
sult, oscillations at the outlet consist of harmonics that are gathered around 
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Figure 11. Characteristic exponents pi and p2 vs. uj. Curve P2{^) has two positive 
maxima indicating the instability of corresponding FDS. Positions of the maxima 
coincide with the frequencies of Hopf and Turing modes that are indicated by the 
dashed lines "H" and "T", respectively, a = 15, 6 = 5.5, 5 = 10, r = 1, = 13. 

maxima of P2{'^)i i-e. around linear modes of the initial system ([1]). This is 
remarkably coincides with our qualitative observations in the above sections: 
when an FDS is unstable inlet noise stimulates the growth of an unstable 
linear mode to an oscillatory solution. 

When the flow rate grows the structure of p2(i^) varies in two different ways. 
The first one is illustrated in Fig. [121 This corresponds to the first scenario of 
the FDS stabilization that is represented in Fig. [TOlfa.b and c) and has been 
discussed in Sec. HI Different curves in Fig. [12] are plotted for different flow 
rates: 

(1) The curve p2 has two positive maxima. Selected noise harmonics excites 
oscillatory mode which destroy the FDS. Hence the is linearly unstable. 

(2) The first maximum moves towards the origin while the second one be- 
comes negative. The Turing mode can not be excited any more, but the 
other one (i.e., the Hopf mode in this figure, but also it can be a DIFI 
mode) destroys the pattern. The system is still unstable. 

(3) The first maximum continues the moving to origin. Thus, frequencies of 
selected noise harmonics differ more and more from the Hopf frequency. 
Because of this detuning the excitation of Hopf oscillations becomes less 
effective. 

(4) Two symmetric maxima (the left one is not shown) merge with a central 
minimum and these three extremums transform into a single maximum 
at the origin. Now all the noise harmonics have non-positive amplification 
ratio and those which are close to the Hopf mode are even suppressed. As 
a result, the noise can not excite oscillations. The FDS is linearly stable. 
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Figure 12. First scenario of the FDS stabilization. Bullet points show the successive 
positions of the right maximum (the left symmetric maximum is not shown) as the 
flow rate grows. Line "H" indicates the Hopf frequency and "T" mark the Turing 
frequencies for different flow rates, a = 15, 6 = 4, 5 = 10, r = 1, (/>! = 7, (1)2 = 8, 
<^3 = 9,(1)4 = 11. 

Because the second couple of maxima passes zero first, it does not play an 
important role in the stabilization. Thus, we classify as the first scenario all 
cases when these maxima absent at all. For example, these cases are shown in 
Fig.[IO]:a,b). 

Fig. [T^ corresponds to the second scenario that has been shown in Fig. [TUT d): 

(1) There are two positive maxima and the FDS is linearly unstable. 

(2) The first maxima merge (only the right part of the symmetric curve is 
shown) but the second ones are still positive and thus responsible for the 
instability of FDS. 

(3) The second couple of maxima pass zero (only the right one is shown in 
the figure). As a result, no oscillations are excited and the FDS becomes 
stable. 

So we see that two scenarios of the FDS stabilization that have been revealed 
from numerical simulations now are confirmed by a linear stability analysis. 

We can provide a direct verification of linear nature of the instability of FDS 
near the critical point 0fds- Let us suppose that the inlet perturbation is nearly 
equal to the saturated FDS amplitude, so that we can neglect the boundary 
layer of an FDS pattern. In the presence of inlet noise with an amplitude 6, 
the FDS occupies an area of a finite size, say h. At the right edge of this area 
the noise reaches the amplitude 6i which is sufficient to destroy the FDS. If 
the instability is Unear, we can write 9i = Oe^^'^ , where p2 = P2(i^max) is the 
characteristic exponent responsible for the stability of FDS. Hence, the size of 
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Figure 13. Second scenario of the linear stabilization of FDS. Bullet points show 
how the second maximum becomes negative as the flow rate grows, a = 20, b = 9, 
S = 20, r = 1, 4>i = 23, (/>2 = 30, 03 = 35. 
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Figure 14. Size of an area of a locahzed FDS h vs. logarithm of the noise amplitude 
6. Observe linear dependance as predicts Eq. (fT2]l a = 20, 6 = 3.5, S = 3, r = 1. 

an FDS area varies linearly with the logarithm of the noise amplitude as 

/i= (In0i-ln^)/p2. (12) 

This equation can easily be verified directly. Solving numerically Eqs. ([1]) 
for different noise amplitudes, we can find corresponding sizes of the FDS 
area. Fig. [H] displays the resulting curves obtained for different flow rates. As 
expected, h decays linearly with In 6'. 

Eq. (fT2l) can be used for experimental measuring of the critical flow rate of 
linear stabilization 0o- A factor in Eq. (JT2i) is the reciprocal characteristic ex- 
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ponent p2(^^max)- To find this factor one needs to vary amplitude of the inlet 
fluctuations and measure the size h of the FDS area. Unfortunately, it is not 
easy to perform this practically. But because we deal with a linear problem, 
each inlet harmonic is amplified separately. So we can substitute the fluctu- 
ations with the most amplified harmonid ^ I. The size h of the corresponding 
FDS area should be equal to that producing by a wide-band inlet noise. The 
frequency of this harmonic is ciJmax- This value may be calculated using the 
Floquet theorem as we have done above for our system. Thus varying the 
amplitude 9 of the inlet forcing and measuring corresponding sizes h one can 
obtain a line h{\n.9) and then calculate p2 as an inverse slope of this line. 
After that one needs to repeat the measuring at different flow rates. Knowing 
P2 vs. flow rate, one can extrapolate the dependance p2(0) to a point where 
P2 vanishes and find the critical flow rate 0o- 



6 Summary and results 

We studied the stability of FDS patterns in presence of a noise at the inlet. 
The following picture was revealed. If the flow rate is not so high, fluctuations 
at the inlet stimulate the growth of an oscillatory mode and thus prevents the 
formation of FDS. For higher flow rates, however, an FDS pattern becomes 
linearly stable. We found two scenarios of the linear stabilization: 1) detuning 
of the most amplified noise harmonics from the resonance with an oscillatory 
mode; 2) stopping of a noise amplification. 

In our previous paper [26] we showed that linearly stable FDS pattern still 
can be destructed by a noise if its amplitude is sufficiently high. A threshold 
value of the noise amplitude depends on the flow rate as a power law. Or 
vice versa: for a fixed noise amplitude there is a flow rate above which an FDS 
pattern becomes stable. This was referred to as a nonlinear instability of FDS. 
In present paper we consider all possible FDSs: 1) FDS is a stationary spatial 
distribution of Hopf oscillations; 2) FDS emerges because of presence a DIFI; 
3) FDS based on the Hopf mode appears when also a Turing mode is present 
and both Hopf and Turing modes have high growth rates; 4) FDS appears 
when the Turing mode is much stronger then the Hopf mode. These numerical 
experiments confirm our previous analysis. Also, the numerical simulations 
reveal two scenarios of the stabilization that remarkably correspond to those 
found within a linear approach. 

Let us consider an equation for a perturbation to FDS. (Now we mean the 



Experiments where FDS is forced by a periodic oscillations at the inlet are re- 
ported by Kffiin and Menzinger [415] and also similar problem is considered by 
Taylor at al. [17j . 
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nonlinear equation whose linearized version has been studied in Sec. |5l see 
Eq. ([8]).) Our observations can be summarized in terms of bifurcations of 
the perturbation equation. A qualitative picture is following. When an FDS 
is linearly unstable, the perturbation equation has a kind of unstable fixed 
point. A bifurcation takes place when the flow rate grows and our analysis 
predicts two scenarios for this. The result of the bifurcation is the split of the 
unstable fixed point into a stable fixed point and unstable limit cycle. This 
corresponds to a nonlinear instability of the FDS. Small inlet noise decays 
because the system remains close to the stable fixed point. But higher noise 
brings out the system above the unstable limit cycle. An amplification of noise 
results in the destruction of FDS. If we increase the flow rate while keeping 
the noise amplitude unaltered, a position of the limit cycle changes and the 
system returns back to a vicinity of the stable fixed point. The unstable limit 
cycle inherits properties of the unstable fixed point, so that two scenarios of 
stabilization are observed when the system passes this cycle. These are the 
scenarios that are obtained in numerical experiments in Sees. [3] and HI 

The destruction of FDS in presence of an inlet noise takes place because this 
periodic pattern serves as a selective amplifier of the forcing. The highes am- 
plification takes places at frequency equal to an unstable oscillatory mode. 
This resonant interaction excites the oscillatory solution that suppresses the 
FDS. Hence, the specific type of an inlet forcing is not very important. The 
sufficient requirement is that its spectrum contain harmonics close the reso- 
nance. In particular, it means that our results remains valid for a Gaussian 
noise. We tried this type of noise. The results were qualitatively similar. 

A pure periodic forcing at the inlet is already studied well. As reported 
in [^lISHTTj . the inlet oscillations of concentrations stimulate the destruction 
of FDS provided that the frequency is close to the natural frequency of a 
system. In the other words, the destruction takes place when the forcing is 
in the resonance with a mode responsible for the formation of FDS. A minor 
decreasing of the inlet frequency from the resonance results in the pulsating 
waves as in [1] or in the zigzag patterns [T7]. These structures are very sim- 
ilar to our intermediate situations represented in Figs. Mjo) and[7](b). These 
results, in fact, corresponds to our first scenario of the FDS stabilization via 
detuning from the resonance. 

Also the resonance with the natural frequency of a system is reported by 
McGraw and Menzinger [15]. They apply a periodic modulation to the fiow 
rate and observe a localized FDS area suppressed downstream by an oscillatory 
solution. It is interesting that the suppression takes place not only at the 
natural frequency. The resonance and further suppression is also observed at 
nonlinear rations 2 : 1 and 3:1! 

Finally, we suggest an idea for an experimental measuring of a fiow rate of 
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linear stabilization. Applying to the FDS inlet a weak periodic forcing with 
specific frequency one can find a characteristic exponent that is responsible 
for the linear instability. One needs to find this exponent for different flow 
rates and then extrapolate the dependance to a point where the exponent 
vanishes. A corresponding flow rate is the sought critical value of the linear 
stabilization. 

This work was partially supported by grant of RFBR (no. 04-02-04011). 
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